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Abstract —We consider a full-duplex decode-and-forward sys¬ 
tem, where the wirelessly powered relay employs the time¬ 
switching protocol to receive power from the source and then 
transmit information to the destination. It is assumed that the 
relay node is equipped with two sets of antennas to enable full- 
duplex communications. Three different interference mitigation 
schemes are studied, namely, 1) optimal 2) zero-forcing and 
3) maximum ratio combining/maximum ratio transmission. We 
develop new outage probability expressions to investigate delay- 
constrained transmission throughput of these schemes. Our anal¬ 
ysis show interesting performance comparisons of the considered 
precoding schemes for different system and link parameters. 

I. Introduction 

Most wireless radios so far have adopted half-duplex (HD) 
communications due to the challenge of handling loopback 
interference (LI) generated from simultaneous transmit/receive 
operation. However, thanks to the progress made on LI 
suppression recently, full-duplex (FD) communications have 
emerged as a viable option El-0- In theory, FD operation 
can double the HD capacity, hence is a key enabling technique 
for 5G systems. 

On the other hand energy harvesting communications is a 
new paradigm that can power wireless devices by scavenging 
energy from external resources such as solar, wind, ambient 
RF power etc. Energy harvesting from such sources are not 
without challenges due the unpredictable nature of these 
energy sources. To this end, wireless energy transfer has been 
proposed as a promising technique for a variety of wireless 
networking applications a. 

RF signals can carry both information and energy and 
pioneering contributions quantifying this fundamental tradeoff 
have been reported. In order to remedy practical issues (same 
signal can not be used for both decoding and rectifying) 
associated with simultaneous information and energy transfer, 
two practical techniques, i.e., time-switching (TS) and power¬ 
splitting (PS) were proposed in O- Both TS and PS apply in 
different network topologies and integration of the RF energy 
transfer into cooperative relay networks is an interesting 
research topic. Different HD relay networks have been studied 
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considering amplify-and-forward (AF) relaying a, a, relay 
selection 0 and multiple antenna relay systems J8]. In 0, 
the achievable throughput of FD AF and DF relaying systems 
with TS has been studied. However, 0 only assumed single 
transmit/receive antenna at the relay. 

Inspired by the FD approach, in this paper we consider a 
source-relay-destination scenario where the multiple antenna 
FD relay is powered via wireless energy transfer from the 
source. The reason for the adoption of multiple antennas at the 
relay is two-fold: (1) employment of an antenna array helps 
the relay to accumulate more energy (2) spatial LI cancellation 
techniques can be deployed. Specifically, we investigate the 
outage probability and the delay constrained throughput by 
considering several precoding schemes at the relay. Our results 
are general in the sense that we consider arbitrary number of 
receive/transmit antennas at the FD relay input/output. 

In summary, the contributions of this work are as follows: 

1) Assuming different precoding schemes at the relay, 
namely, optimal, zero-forcing (ZF) and maximum ratio 
combining (MRC)/maximum ratio transmission (MRT), 
we develop new expressions for the system’s outage 
probability, which are helpful to investigate the effect of 
key system parameters on performance metrics such as 
the outage probability and delay-constrained throughput. 

2) In the case of ZF, we present simple high signal-to- 
noise ratio (SNR) expressions for the outage probability 
which enable the characterization of the system’s diver¬ 
sity order and array gain. Moreover, we compare the 
performance of FD and HD modes to show the benefits 
of FD operation. 

II. System Model 

We consider a DF relaying network consisting of one source 
S, one relay R, and one destination D. Both S and I) are 
equipped with a single antenna, while R is equipped with Mr 
receive (input) antennas and Mr transmit (output) antennas to 
enable FD operation. We assume that the S to D link does 
not exist, due to severe shadowing and path loss effect. 

It is also assumed that the relay has no external power 
supply, and is powered through wireless energy transfer from 
the source. We adopt the time-sharing protocol 0, hence 
the entire communication process is divided into two phases, 
i.e., for a transmission block time T, a fraction of the block 
time is devoted for energy harvesting and the remaining time. 



(1 — a)T, is used for information transmission. It is also 
assumed that the channels experience Rayleigh fading and re¬ 
main constant over the block time T and varies independently 
and identically from one block to the other. 

During the energy harvesting phase, the received signal r e 
at the relay can be expres sed as 

( 1 ) 

where Pr is the source transmit power, d\ is the distance 
between the source and relay, r is the path loss exponent, 
hg/{ is the Mr x 1 channel vector for the S-R link, i.e., 
input antennas at R are connected to the rectifying antenna 
(rectenna), x e is the energy symbol with unit energy, and n r 
is the zero mean additive white Gaussian noise (AWGN) with 
unit variance. We assume that the energy collected during the 
first phase is fully consumed by the relay to forward the source 
signal to the destination. Hence, the relay transmit power can 
be computed as P r = ^Pg||hg/;|| 2 , where n = with p 
denoting the energy conversion efficiency. 

Now, let us consider the information transmission phase. 
The received sig nal a t R can be written as 

[Ps 

r[n] = \ — h srx s [n ] + H RR x R [n] + n R [n] , (2) 

y “i 

where xs[n\ is the source information symbol with unit 
energy, and x#[n] is the transmitted relay signal satisfying 
E |xfl[n]x^[n] | = P r , and H## denotes the Mr x Mr LI 
channel. Upon receiving the signal, R first applies a 1 x Mr 
linear combining vector w r on r[n] to obtain an estimate 
of xs, then forwards signal to the destination D using the 
Mt x 1 transmit beamforming vector w t . It is assumed that 
IKII = ll w r|| = L 

The relay 's e stimate xs[n] = w r r[n] can be expressed as 

fP ,s 


x s M = 1 / -Jf W> hgflXg [n\ 
a 1 


-w r H flfl x fl [n]+w r n R [n]. (3) 


end SINR in ([5]) is maximized. Specifically, for a fixed value 
of a, the SINR maximization problem can be formulated as 
max 7 (in Eq. ©) ( 6 ) 


W t ,W r 

S.t. 


= l. 


II W HI = ||w t 

In order to solve the problem in ©, we first fix w t and 
optimize w r to maximize 7. Therefore, the optimization 
problem can be re-formulated as 

|w r hs fl | 2 


max 


s.t. 


ItPs 

dT 


||hsit|| 2 |w r H M w t | 2 

11 = 1 , 
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(7) 


which is a generalized Rayleigh ratio problem. It is well known 
that the objective function in {7} is globally maximized when 
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( 8 ) 


Accordingly, by substituting w r into the objective function, 
the optimization problem in © can be re-expressed as 


. , l.. u ||2 ^||hg«|| 2 ||h^H M w t || 2 

max Fs mm — hs#- - - ^ -■ 

l dT 11 " l + ^||h Si? || 2 ||H^w t || 2 ^ 


kPs 


did 


-||h S ij|| 2 |hijz3W t |' 


s.t. 


||w t || = 1, (9) 

which is still difficult to solve. Therefore, instead of © 
we solve the following problem by introducing an auxiliary 
variable t, as 

|h flD w t | 2 (10) 


/(f) = max 

II w t ||=i 


s.t. 


K-Pfl 


£||hg*|| 2 ||h^H^w t || 2 
l + ^IM 2 ||H Kfl w t || 2 


= f. 


The relay transmit signal is given by x^[n] = 

y/P^WtXs [n — <5] where <5 accounts for the time delay caused 
by relay processing. Finally, the received signal at D is 

Vd M = \ -Tf x R [n\ + n D [n\ . (4) 

V ^2 

where lino is the 1 x Mt channel vector of the 11 I) link, 
d 2 is the distance between the relay and destination, and iid 
is the zero mean AWGN with unit variance. 

With the DF protocol, end-to-end signal-to-interference- 
plus-noise ratio (SINR) can be written as 

7= (5) 

P s |w r h s _R | 2 kP s .. ||2| |2 \ 

5-77 T -- 7 J- 12. Jr '^ FJF \\ h SR \\ hflflWt . 

III. Joint Precoding/Decoding Designs 

In this section, we consider several precoder/decoding de¬ 
signs to suppress/cancel the effect of LI at the relay, each of 
which offers different performance-complexity tradeoff. 

A. The Optimal Scheme 

In this subsection, our main objective is to jointly design the 
precoder and the decoder at the FD relay so that the end-to- 


This is a nonconvex quadratic optimization problem with 
quadratic equality constraint. To solve the problem in ( ITOl . we 
apply a similar approach as in ED to convert the optimization 
problem to 

max tr(hJj£,Wfh/?x>) (H) 


w t 


s.t. tr(W 1 /H^ fl hsnh]; i jH R . R 

tr(Wt) = 1, 

where W t = w t wj is a symmetric, positive semi-definite 
matrix. In order to solve (DU. we can resort to the widely 
used semidefinite relaxation approach. By dropping the rank- 
1 constraint, the resulting problem becomes a semidefinite 
program, whose solution W t can be found by using the 
method provided in ltl3l Appendix B]. 

Denoting the optimal objective value of (fTTb as /(f), the 
SINR maximization problem can be formulated as 

max min (p s (Ll|h s rf-() ,|L||h SB ||V(f)) . (12) 

Therefore, in order to solve ©, it remains to perform a one¬ 
dimensional optimization with respect to the variable f. 
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B. Transmit ZF (TZF) Scheme 

In the transmit ZF scheme, the FD relay takes advantage of 
the multiple transmit antennas to completely cancel the LI. To 
ensure this is feasible, the number of the transmit antennas at 
R should be greater than one, i.e., Mt > 1. In this case w r = 
ira and Wf is the solution of the following optimization 
problem: 

max |h_R£)W t | 2 (13) 

W t 

s.t. ||w t || = 1 , hl R U RR w t = 0 . 

From the ZF constraint, we know that w t lies in the null space 

of hl R H RR . Denoting B = I - we have 

on II^srHrrII 2 


where Q(a, x) = T(a,x)/T(a) and T(a,x) is the upper 
incomplete Gamma function EH Eq. (8.350.2)], and Pl = P s 
is the SNR of the first hop. 

To the best of the authors’s knowledge, the integral in (IThl) 
does not admit a closed-form expression. However, this single 
integral expression can be efficiently evaluated numerically 
using software such as Matlab or Mathematica. 

To gain further insights, we now look into the high SNR 
regime and derive a simple approximation for the outage prob¬ 
ability, which enables the characterization of the achievable 
diversity order of the TZF scheme. 

Proposition 1. In the high SNR regime, i.e., p\ —> oo, the 
outage probability of FD relaying with the TZF scheme can 
be approximated as Ell at the top of the next page, where 
ip(x) is the digamma function HlO\ Eq. (6.3.1)]. 


C. Receive ZF (RZF) Scheme 

As an alternative solution, the transmit beamforming vector 

h* 

can be set using the MRT principle, i.e., w t = |i h R - P |i , and 
w r based on the ZF criterion. To ensure feasibility of ZF, R 
should equipped with M R > 1 receive antennas. 

By using similar procedure as shown for the transmit ZF 
scheme, the combining vector w r can be obtained as w^ F = 


Dh sr 

IlDhsfiH 


with D 


IHrrIirb) || 2 


rtn 


D. MRC/MRT Scheme 


For the MRC/MRT scheme, w r and w t are set to match 
the first hop and second hop channel, respectively. Hence, 


r MRC 


= & and = 


||h KD || ' 


IV. Outage Probability 

In this section, we investigate the outage probability of 
the considered FD relay system assuming TZF, RZF, and 
MRC/MRT schemes. In case of the optimal scheme, derivation 
of the outage is difficult and we use simulations in Section [V] 
The outage probability is an important performance metric, 
which is defined as the probability that the instantaneous SINR 
falls below a predefined threshold, 7 th. Mathematically, it can 
be written as 

Pout = Pr(7 < 7th) = -F 7 (7th)- (14) 


A. TZF Scheme 

By substituting the xvf F and w^ RC into <0, the end-to-end 
SNR 7 tzf can be expressed as 

P S lit m 2 , 


7TZF = 


(1 — a)d[ 




-Jr\\h R D \\ 2 


(15) 


where h RR > is an ( Mt — 1) x 1 vector with ||/i,r.d|| 2 ~ 
* 2 (m t — l)- Let Y = min ( x - If II^hdII 2 )- Then invoking 
the cumulative density function (cdf) of Y presented in [?, 
Appendix II], the cdf of 7 tzf can be obtained as 


^TtzfO) = 1 - 


M R —1 -X 


£1 1 T(Mr) 


e X dx. (16) 


Proof: Due to limited space proof is omitted. ■ 

By inspecting ( 1 17 l i. we see that the TZF scheme achieves a 
diversity order of min (M R , M R — 1). Moreover, we notice that 
for M r + 1 = Mt, decays as pf MR ln(pi) rather 

than Pi Mr as in the conventional case which implies that in 
the energy harvesting case the slope of F 7 ~ f (^) converges 
much slower compared with that in the constant power case. 


B. RZF Scheme 

Invoking 0, and using w f f and w f MRT , the end-to-end SNR 
7rzf can be expressed as 

7 rzf = ZPs min —||/iad|| 2 ^ , (18) 

where Z ± ± (||/rsfl|| 2 + l^i| 2 ) and X x = 

It is well known that Z follows central chi-square distribution 

with 2 M R degrees-of-freedom, denoted as Z ~ xi.M,; an d that 

X\ follows a beta distribution with shape parameters M R — 1 

and 1, denoted as Xi ~ Beta(M R — 1,1) with Fx 1 {x) = 

x M r ~\ 0< x < 1, E2 P- 138], 

Moreover, by denoting Y\ = -Tr\\h R T,\\ 2 , we have FyAy) = 

a 2 

P (Mt, where P(a,x) = 'y(a,x)/T(a) is the regu¬ 

larized lower incomplete Gamma function EOl Eq. (6.5.1)]. 
With F Xl (x) and Fy 1 ( y ) in hand, the cdf of 7rzf can be 
expressed as ( | 1 9| i. Although Eq. ( 1 1 9| ) does not admit a closed- 
form solution, it can be efficiently evaluated. Now, we look 
into the high SNR regime, and investigate the diversity order 
achieved by this scheme. 


Proposition 2. In the high SNR region, i.e., p\ —> oo, the 
outage probability of the FD relaying system with the RZF 
scheme can be approximated as 




r (Mr) 


T(Mr) 


(31) 


Mr —1 


1 + 


M r < Mt +1, 

r(M T ) J \ pi ) 

Mr = Mt + 1 , 


T(Mr-Mt) ( dl\ Mt ( d(z\ 
, r(MR)r(M T +i)\ k J \ pi J 


Mt 


M] j > Mt + 1 . 


Proof: The proof is omitted due to space limit. ■ 

Proposition [2] indicates that the RZF scheme achieves a 
diversity order of min(MR — 1 ,Mt). 
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i_ | i v C-D fc+1 (dl \ MT+k ~ 1 _ 

r(Mfi+l) ^ r(M r -l)r(M R ) ^ k\(k+M T ) \ K ) m r - 

^Mr-Mt+1 

r(MH+i) + TTO ( ln (^i) - ln + V>(1)) (£) *) (tt) *> 

r(M H -M r + l) /d ££\ MT_1 

r(M T )r(M B ) \ ^ \pi) ’ 




Mt > Mr +1, 

Mr = Mr +1 
Mr <C AT/j -|- 1 


(17) 


R mzF (z)=l-Q(Af« 


pi 
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+ r(M fl ) 




d[d££\ 

Kpi X/ 




d[d£A 

K/?l X/ 



(19) 


C. The MRC/MRT Scheme 


Substituting w^ RC and w^ RT into 0. analysis of the 
MRC/MRT scheme for arbitrary Mt and Mu appears to be 
cumbersome. Therefore, in the sequel we consider two special 
cases as follows: 

I ^srHrrH^p I . 


Case-1) Mt = 1, Mu > 1: In this case 
given by 


IlhsRhflD || 2 


X! = K 


-h 


hsR 


rr\ = 




\h 


ititll 


( 20 ) 


For notational convenience, we define c\ = -pf, C 2 = Psk ° t rr , 
C 3 = -i-rff and write the end-to-end SINR as 


7m RC =mm 


Clllhsitll 2 


^IhsRT^WhRRW 2 + 1 


'C3\\hsR\\ 2 \hRD\' 


Let us denote X = C\/ (^ 2 X 1 + and Y = C 3 Y 1 Y 2 , with 

Y\ = \\hsR\\ 2 and Y 2 = |/i_r_d| 2 . Conditioned on Y\, the 
random variables (RVs) X and Y are independent and hence 


=1 - r f Xi - ~ 

Jjl- \C 2 \Z y 


x 1-4 


fvi(y)dy. (21) 


(£). 

In order to evaluate (EH we require the cdf of the RV, X-\. 
Note that ||/i fl /{|| 2 ~ xIm R ’ Z ~ ■p^pr is distributed as 
Z ~ Beta(l, Mr — 1), and the cdf of Xi can be evaluated as 

1,Mr ^ ( 22 ) 


where G™ 

(9.301)]. Now, using the cdf of RV K>, and substituting 
into (EH we obtain 


(z- I “ 1 ;;;“^ is the Meijer G-function JTT] Eq. 


^7MRC (~) 1 


1 


G. 


21 


T{M r )J a ~ 23 
C 1 

X y MR ~ 1 e~^ V+T rs) dy. 


1 ,Mr 

Ij Mr, 0 


(23) 


To the best of the authors’ knowledge, the integral in (EH 
does not admit a closed-form solution. However, (EH can be 
evaluated numerically. 

Case-2) Mr = 1, M t > 1: In this case 


simplifies to 

Yi = (hflfl*tdiag{l,0 ) --- ,0}*^) = \h RRA \ 2 , (24) 


where h,RR = tiRR&f Hence, 7 mrc can be written as 
^ cxlhsRl 2 


7m rc = mm 


C2 | h,SR | 2 \tlRR,i | 2 + 1 


,C3\hsR\ || ft/to || ^ 


Let us define X = C\Xi/ (c 2 X\Y\ + 1), and Y = C 3 X 1 Y 2 
where X\ = \1isr\ 2 , V 2 = ||/i/to|| 2 - Note that conditioned on 
Xi, the RVs, X and Y are independent and hence we have 

/»oo 

FW(~) = i-y (l-^xix^z)) (l-F YlXl (z)) f Xl {x)dx, 

= 1 -J (l ^M T ,Mje~ x dx. (25) 


D. Half-Duplex Scheme 

We now present the outage probability of the HD relaying 
system, which serves as a benchmark for performance com¬ 
parison. The energy harvesting phase of the HD relaying is 
the same as that of the FD relaying system. However, in the 
information transmission phase, the remaining (1 — a) portion 
of block time is equally partitioned into two time slots for 
source and relay transmissions 0. The end-to-end SNR of 
the HD relaying scheme can be computed as 

7hd = ||ft.s/t|| 2 min(ci,2c 3 ||/i fl D|| 2 ) ■ (26) 


The required cdf can be obtained by replacing Mt — 1 with 
Mt in (IT6l) . and is given by 


-^7hd ( z ) — 1 


/; 




d l z T{Mr) 


X M R -l e ~ x dx _ 


(27) 


V. Numerical and Simulation Results 
We now present numerical and simulation results to investi¬ 
gate the impact of key system parameters on the performance. 
In all cases, we have set rj = 1. 

Rg. m compares the outage probability of the considered 
schemes with different antenna configurations and for a = 0.5, 
di = g ?2 = 1, < 7 r R = 0.1, and 7 th = 0 dB. When the transmit 
power is high and a remains fixed, an excessive amount of 
energy will be collected at the relay, which is detrimental for 
the MRC/MRT scheme since it results in strong LI. Therefore, 



























Fig. 1. Outage probability versus pi of the TZF, RZF and MRC/MRT 
schemes for different antenna congurations. ( 7 ^ = 0 dB, and a = 0.5). 



Fig. 2. Delay-constrained throughput of the proposed precoding schemes 
(Mr = Mr = 4, pi = 10 dB. &\r = 0.3, d\ = d ,2 = 2, r = 3.1 and 
R c = 1 bit/sec/Hz). Results for the MRC/MRT scheme are from simulations. 

the outage performance of the MRC/MRT scheme exhibits 
an outage floor at high SNRs. On the contrary, the outage 
probability of the TZF and RZF precoding schemes decays 
proportional to the diversity orders provided in Proposition!]] 
and [2] respectively since the relay is capable of canceling LI. 
However, a proper choice of a can improve the outage of the 
MRC/MRT scheme to some level at high SNR. On the other 
hand observation of different curves in the low SNR region 
reveal that the MRC/MRT scheme outperforms the TZF and 
RZF schemes even at high LI level. This is because at low 
SNR, overall interference can be treated as noise and therefore 
MRC filtering helps to maximize the SNR. Comparing the 
TZF and RZF schemes with the same diversity orders and 
different receive antenna numbers (i.e., TZF, with Mr = 2 
and Mt = 2, and RZF with Mr = 3 and Mr = 1) we see 
that the additional receive antenna could harvest more energy 
to facilitate information transfer. Moreover, for the case where 
Mr = Mr, TZF achieves a higher array gain. 

Fig. U shows the impact of optimal time split a on the 
delay-constrained transmission throughput defined as: R(a) = 
9(1 — P out )f? c (l — a), where 6 = 1 for FD and 9 = 0.5 for 
HD and the optimal time portion a can be obtained from 
a* = argma xR(a), 0 < a < 1. 


As expected, the optimal scheme exhibits the best throughput 
out of all precoding schemes studied. The superior perfor¬ 
mance of the optimal scheme is more pronounced especially 
between 0.4 and 0.8 values of a. The highest throughput 
with optimized a for the optimal, RZF, MRC/MRT and 
TZF schemes are given by 0.382, 0.374, 0.358 and 0.315, 
respectively. Moreover, we see that each one of TZF, RZF and 
MRC/MRT precoder designs can surpass the others depending 
on the value of a. Finally, all schemes achieve significant 
throughput gains as compared to the HD mode. 

VI. Conclusion 

In this paper, we have studied the outage probability and 
throughput of FD MIMO relaying with RF energy transfer. We 
designed the optimal precoder/decoder as well as investigated 
several sub-optimal low complexity precoding schemes. The 
MRC/MRT scheme provides a better outage performance at 
low-to-medium SNR, while the ZF precoders outperform the 
former at high SNR. Further, it is demonstrated that all 
proposed FD precoding schemes attain significant throughput 
gains compared to the HD mode. Therefore, FD relaying 
is a promising solution for implementing future RF energy 
harvesting cooperative communication systems. 
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